Quantum Superposition State Production by Continuous Observations and Feedback 
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We present a protocol for generation of superpositions of states with distinguishable field ampli- 
tudes in an optical cavity by quantum non-demolition photon number measurements and coherent 
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By suitably tailored optical pulses it is possible to co- 
herently manipulate the states of small quantum systems 
and for example to steer molecular processes and chem- 
ical reactions. Methods and concepts from this research 
have spread to the field of quantum information theory 
which, even with quantum error correction, requires a 
very high degree of control [l|. As an example, quan- 
tum optimal control techniques can substantially improve 
the performance of elementary quantum gates with cold 
neutral atoms Q. Optimal control methods aim at ma- 
nipulating a few external parameters, e.g., currents and 
magnetic fields of an atomic trapping potential, in such 
a way that an initial state of the system evolves into 
the desired final state with high fidelity. These tech- 
niques are open- loop, i.e., they do not exploit the knowl- 
edge that one can get by observing the system and using 
the measurement outcome in a suitable feedback. Even 
quite simple measurements display powers which are hard 
to match with controllable interactions in terms of the 
states accessible. For example, optical probing of spin- 
polarized macroscopic atomic samples has been used to 
enable atomic spin-squeezing J3[ , entanglement, quantum 
storage Q and teleportation [HJ , and measurements of the 
phase of light transmitted through a modest cavity has 
been proposed as a means to project product states of 
atoms in the cavity into entangled states and to imple- 
ment quantum computation [6|. 

The natural next step is to apply feedback continuously 
in time using the information acquired in real time with 
the measurements. The theory for continuous measure- 
ments and feedback [7|, |8[ combines the non-deterministic 
elements of quantum trajectories with stochastic 

calculus. While these theories describe correctly the out- 
come of a given measurement and feedback scheme, it is 
still an open problem how one identifies reliable schemes 
for a given task. A scheme for photon Fock state genera- 
tion in a cavity has been proposed recently by Geremia in 
Ref. 11] and an analysis of the stability of the feedback 
by Yanagisawa [l2j]. In this Letter we propose a strat- 
egy to generate an equal superposition of two quantum- 
mechanical states with distinguishable field amplitudes 
\S&) oc \A) + | B). Superpositions of states, which are 
well localized in separate regions of the effective position- 



momentum phase space of the field variables, have been 
proposed as useful resources for quantum computation 
[lH and quantum metrology [14]. Other methods for 
generation of such superposition states have been demon- 
strated with light [H, [3] and with trapped ions 
Our principal idea is to apply the method described in 
Ref. pjj to approach a target Fock state, occupying a ring 
in the field amplitude phase space. Before we reach this 
state, the phase space distribution is of crescent shape 
[see Fig. 1(b)], and we feed coherent radiation into the 
cavity to displace the state [Fig. 1(c)]. We subsequently 
start probing the photon number again. In a pictorial 
representation, this scheme will select a quantum state 
with phase space support at the overlap of a new Fock 
state ring and the displaced crescent distribution, i.e., at 
two crossing regions, and hence a linear superposition of 
two quantum states with distinguishable field amplitudes 
may result from the protocol [Fig. 1(d)]. 

A quantum non-demolition measurement of the photon 
number h in a single cavity mode can be accomplished 
by measuring the phase shift of a probe laser field that 
couples via a cross-Kerr effect to the cavity field when it 
passes through an atomic gas inside the cavity. In our 
simulations, we apply physical parameters of a dark- state 
mechanism in the gas for effective coupling of the fields 
taken from Ref. [11], but the formalism is general and 
describes also the effective coupling based on the collec- 
tive atomic motion inside the cavity, studied in Ref. [13] • 
The interaction Hamiltonian is proportional to the prod- 
uct of the cavity and the probe photon number oper- 
ators and it causes a phase shift of the probe without 
exchanging photons with the cavity field. It thus en- 
ables a non-demolition interaction which will gradually 
cause a narrowing of the photon number state distribu- 
tion. Our simulations will apply the stochastic master 
equation technique @, EH, but in order to explain this 
method and to get useful insight, it is worthwhile to es- 
tablish a simple physical picture of the underlying prob- 
ing dynamics. For this purpose, consider the probing 
beam as composed of a succession of segments of dura- 
tion At. The field is in a coherent state, and hence factors 
in a product state of coherent states occupying each seg- 
ment of the beam. The continuous measurement on the 
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probe beam after interaction with the cavity field now 
separates in the detection on each individual segment of 
light, and assuming an incident coherent state with a real 
field amplitude, the phase shift is registered by balanced 
homodyne detection of the phase-quadrature component 
by means of interference with a local oscillator field. 

The interaction between the cavity field and a sin- 
gle segment of light is governed by the unitary operator 
JJ = e~ lM nnp , where M' is a measure of the coupling 
strength. We assume that each segment of the probe 
beam is in a coherent state with a large mean number 
of photons $A£, where $ is the photon flux. There- 
fore, we can write d p — V &At + a' p and expand the 
time evolution operator to lowest order in the quantum 
fluctuations, U = e ~ i(3fl& ^ where (3 = 2 M' V$At : and 
x p = (a'p + a'^)/2. The interaction implies a displace- 
ment of the probe p quadrature proportional to the cav- 
ity photon number, and it is thus useful to expand the 
joint cavity and probe field state in the corresponding 
eigenstate basis, |^ in ) = ^2 n c n \n) J dp p ir~ 1/4 e~ p p /2 \p p ), 
so that the state after interaction becomes 

|*out) = 5>» / dppir-^e-^-^ 2 / 2 ^). (1) 

n 

At this stage, p p is measured, an arbitrary outcome is 
obtained according to the probability distribution, 

f(p P )=^ 1/2 J2\ c n\ 2 e- (p ^ n) \ (2) 

n 

and the projection of (pQ) on the corresponding p p eigen- 
state, yields the updated state of the cavity field, 

\MPp)) = W.ftO E c » e-0»»-fp) a /2 | n> , (3) 

n 

where J\f(/3,p p ) is a normalization factor. When we 
model continuous probing, the parameter /3 is infinites- 
imal, and the effect of the interaction is merely to shift 
the Gaussian probability distribution for p p by /3(n). The 
update of the cavity field state after measurement is also 
infinitesimal due to the weak n dependence of the expo- 
nential factor in (|3|). The random detection can be mod- 
elled by a Wiener noise process, and the conditioned dy- 
namics under continuous measurements can be brought 
on the form of an Ito stochastic master equation (SME) 

dp{t) =MV[h]p(t)dt^^/M^H[h}p(t)dW(t), (4) 

with V[X]p = XpX^ - 1/2(1+ Xp + pX^X) and 
H[X]p = Xp + pX^ - Tr[(X + X*) p] p. In | M = 
2 M' 2 $ denotes the measurement strength and n G [0, 1] 
represents the quantum efficiency of the detection. The 
innovation process, i.e., the difference between the actu- 
ally observed p p and its quantum mechanical expectation 
value with the current quantum state of the cavity field, 



is described by a Wiener process [18], dW(t). This dif- 
ference is due to the shot noise in photo detection. 

While the continuous measurement described above 
will eventually collapse the system on one of the Fock 
states present in the initial state, an important initial 
step of our scheme is to evolve the system towards a 
given Fock state. In Ref. [11], Geremia showed how to 
use the information gradually obtained about n by the 
detection record to feed coherent radiation into the cav- 
ity which in- or decreases the total photon number in 
a controllable manner. This feedback is described by 
adding to (|4]) terms that describes evolution under the 
Hamiltonian H^b (t) = G ef x. Here x = (a + a) ) /2 is the 
cavity field quadrature operator, G is the feedback gain 
factor, and e/ is the feedback policy function, that we 
take to depend on appropriate expectation values for the 
field. A natural choice for the feedback policy function is 
ef((n)) — n* — (n), where n* is the desired photon num- 
ber [ll|]. The feedback Hamiltonian causes a displace- 
ment of the field quadrature operator p = (a — a))/(2i). 
The feedback is proportional to n* — (h) , and a state with 
negative (p) can be shifted to larger negative values of (p) 
and hence typically a larger (n), if desired. 

The functioning of our complete scheme is illustrated 
in Fig. [TJa-e). We find heuristically (for n* c± 10) 
that the following protocol has a high success rate: At 
time t = we start with the electromagnetic vacuum 

(a) and we apply the Geremia Fock state feedback pro- 
tocol towards n*, but only until the quadrature (p) < 
(n* — l) 1 / 4 — yfn*. At that time the Husimi Q-function 
looks like a crescent in position-momentum phase space 

(b) . The state is now shifted towards the positive value 
(p) =0.9 (c) and at that time we start again to 
observe (without feedback) the photon number in the 
system until the desired state is found (d). Finally, the 
state is shifted to a symmetric state around the centre of 
the phase space (e). We stop the action of the probing 
field when the maximum value of the Husimi Q-function 
along the p quadrature axis is below some fixed value S 
(typically 5 < 0.005). 

The protocol does not work in every run of the simula- 
tion/experiment. When we start the measurement on the 
crescent shape state it sometimes happens that the state 
collapses into a quasi coherent state [19| , and sometimes 
the state dynamics becomes unstable. Since we have ac- 
cess to the density matrix conditioned on the measure- 
ment outcome, we know if the collapse takes place, and 
this problem is partially solved by starting over again 
to generate a new crescent state with the number state 
feedback generator. The instability problem is solved 
in our simulations by applying a not too fast ramp of 
increased measurement strength M(tY controlled in the 
experiment by the probe laser power [ll| . In our protocol 
it is important that the probing laser is not switched on 
too fast, but it is also important that it is switched off 
quickly to not perturb the state when it has been created 
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Figure 1: (Color online) Time evolution of the Husimi Q- 
function for a single successful realization of the linear super- 
position state (rj = 1, k = 0, and n* = 10). In (f) the positive 
x component of the Husimi Q-function in (e) is magnified. 
The white ring in (f) indicates a contour for the choice of dis- 
placed squeezed state |a, C) m © that yields the best overlap 
with our state. The lower panel shows the histogram of the 
fidelity for 135 simulations with n* — 10, k = 0, and 77 = 1. 



(Switch-on ~ 200 ns, Switch-off ~ 1 - 2 ns for M = 2.12 
MHz). Using these strategies and allowing a maximal 
production time of 10 M -1 , we observed a success prob- 
ability for the production of a quantum superposition 
state of 51% for n* = 5, 69% for n* = 8, and 77% for 
n* = 10. Those results are obtained in the ideal situation 
of perfect detector efficiency and no cavity decay. 

Our protocol does not automatically favor a superposi- 
tion of coherent states, and we have investigated to which 
extent, the state produced can be written as a superposi- 
tion of two Gaussian, minimum uncertainty states in the 
position- momentum phase space. As a way to quantify 
the quality of the state we use the optimal overlap fidelity 
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max{(*(C,a,0)|p|*(C,a,0)>}, 



(5) 



where |*(C,a,0)) = A/"(C, a, 4>) (|a,C> 
and where ( = re 1 6 is a squeezing parameter, a = x + i y 
is a displacement parameter, and AT(C, ol, (j)) is a normal- 
ization factor. The optimal superposition state parame- 
ters (a, C, (j>) can be determined from the detection record 
in every run of the experiment. The results of 135 at- 
tempts to produce such quantum superposition states are 
summarized in the lower panel of Fig. [TJ which provide 
an average fidelity of about 90%. The success probability 
for n* = 10 does not change appreciably when a cavity 



decay rate of k = 0.005 M and finite detector efficiency 
77 = 0.8 are taken into account in our simulations, but the 
average fidelity decreases to 70%. Note that in our fam- 
ily of test superposition states we allow a relative phase 
(j). This phase is known to the experimenter based upon 
the detection record, but it is not under straightforward 
experimental control. 

We will now explain our quantitative findings. We 
have found numerically that the crescent shaped state of 
Fig. [If b), which is produced with a high success probabil- 
ity, is very close to the so-called crescent state of Ref. [lo[ • 
These states are eigenstates of the non-Hermitian opera- 
tor n — 2 i |£| x [20[. When |£| — > the eigenstate is close 
to a Fock state, while for |£| ^> 1 it is a coherent state. 
We now insert this state in Eq. (pQ) and we can semi- 
analytically follow the effect of probing, without apply- 
ing any feedback. Since the measurement is of quantum 
non-demolition type, the equations (1-3) discussed for 
infinitesimal temporal segments, also apply for the accu- 
mulated effect of measuring on the probe field for an ex- 
tended period of time. The parameter (3 is then larger, so 
that the probability distribution ([2]) for the p p observable 
for a long time interval is no longer well approximated by 
a Gaussian. Instead, p p provides a weak measurement of 
n, which according to Eq. (pQ) would become a projection 
in the limit (3 — » 00. Figure 2(a) shows a plot of the p p 
probability distribution f(p p ) obtained for a typical prob- 
ing time, corresponding to /3 — 0.2, for the crescent state. 
We identify three regions on the curves: (/) a maximum 
for small p p , which is rather independent of n*, (II) a 
central region where f(p p ) decreases, and a second 

maximum at p p ~ n* followed by a rapid decay. Since p p 
measures the cavity photon number n, the appearance 
of the three regions can be qualitatively understood in 
terms of the n distribution of the state or, more illus- 
trative, in terms of the Husimi Q-function. Fock states 
are ring-shaped in phase space and in Fig. 2(c) and 2(d) 
we plot the displaced crescent state together with the 
region borders identified from Fig. 2(a) for n* = 8 and 
for n* = 10. It is seen that region (J) corresponds to 
the part of the state close to the origin. By increasing 
n*, the weight in this region will only change slightly. In 
region (II) the "arms" of the crescent state are more or 
less radial in phase space while in region (///) they come 
together again. When n* is increased, the crescent state 
becomes larger and the border between region (II) and 
(///) must move to larger radii, i.e. larger p p . 

From the phase space plots it is to be expected that 
region (II) is where a weak n measurement will cut out 
two well separated portions of the Husimi Q-function. 
This is confirmed by Fig. 2(b), which shows the maximum 
value of the post-measurement Husimi Q-function on the 
p-axis as a function of the measured p p . We see that in 
the central region (J/), the state produced has very small 
values of the Husimi Q-function on the p-axis, while the 
surrounding regions (/) and (III) are clearly not useful 
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Figure 2: (Color online) (a) Probability distribution f(p P ) 
for the outcome of the measurement of the probing field 
quadrature p p . (b) Maximum value of the Husimi Q-function 
Q(x = 0,p) plotted as a function of the measurement outcome 
Pp. The black (continuous) lines correspond to n* = 8 while 
the red (dashed) ones to n* — 10. In (c) and (d) are shown 
the corresponding Husimi Q-functions of the displaced cres- 
cent state together with the region boarders identified from 
(a) and (b). In all pictures f3 = 0.2 has been taken. 



for our purpose. Since the second maximum in Fig. 2(a) 
moves further to the right with increasing n*, the central 
region (II) with the successful outcomes becomes larger 
and this explains that our protocol works better for larger 

It should be noticed that the semi-analytical model 
based on the crescent state [lo| suggests that the opti- 
mal superposition state parameters (a, £, (j>) be uniquely 
determined by n*, the duration of the final probing stage, 
and the corresponding integrated signal p p . In an experi- 
ment, this is much easier than the numerical optimization 
based upon the solution of the SME (|4|). 

In conclusion we have proposed to generate non- 
classical states of light in a cavity by using quantum 
measurements and feedback. A protocol for production 
of highly non-classical superposition states with high fi- 



delity and success probability was proposed. We are cur- 
rently working on the generalization of the ideas pre- 
sented in this Letter to the generation of similar states 
of atomic ensembles. 
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